Density estimation by total variation penalized likelihood driven by the sparsity `1 information criterion

We propose a non-linear density estimator, which is locally adaptive, like wavelet estimators, and positive everywhere, without a logor root-transform. This estimator is based on maximizing a non-parametric log-likelihood function regularized by a total variation penalty. The smoothness is driven by a single penalty parameter, and to avoid cross-validation, we derive an information criterion based on the idea of universal penalty. The penalized log-likelihood maximization is reformulated as an l1-penalized strictly convex programme whose unique solution is the density estimate. A Newton-type method cannot be applied to calculate the estimate because the l1-penalty is non-differentiable. Instead, we use a dual block coordinate relaxation method that exploits the problem structure. By comparing with kernel, spline and taut string estimators on a Monte Carlo simulation, and by investigating the sensitivity to ties on two real data sets, we observe that the new estimator achieves good L1 and L2 risk for densities with sharp features, and behaves well with ties. SARDY, Sylvain, TSENG, Paul. Density Estimation by Total Variation Penalized Likelihood Driven by the Sparsity l1 Information Criterion. Scandinavian Journal of Statistics, 2010, vol. 37, no. 2, p. 321 337 DOI : 10.1111/j.1467-9469.2009.00672.x